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Abstract. We use exact numerical diagonalization in the search of fractional quantum Hall 
states with non-Abelian quasiparticle statistics. For the (most promising) states in a partially 
filled second Landau level, the search is narrowed to the range of filling factors 7/3 < v e < 8/3. 
In this range, the analysis of energy spectra and correlation functions, calculated including finite 
width and Landau level mixing, supports the prominent non-Abelian candidates at v e — 5/2 
(paired Moore-Read "pfafian" state) and 12/5 (clustered Read-Rezayi "parafermion" state). 
Outside of this range, the noninteracting composite fermion model with four attached flux 
quanta is validated, yielding the family of quantum liquids with fractional, but Abelian statistics. 
The borderline v £ = 7/3 state is shown to be adiabatically connected to the Laughlin liquid, 
but its short-range correlations are significantly different. 



1. Introduction 

The incompressible quantum liquids (IQLs) p] continue to be the subject of extensive studies 
ever since the famous discovery of the fractional quantum Hall (FQH) effect |2j. The IQLs 
are formed by two-dimensional electrons placed in a high magnetic field B which causes them 
to fill a particular fraction v of one of the lowest Landau levels (LL n , n = 0, 1, ...). The 
most recent storm of interest in the IQLs is motivated by the concept of "topological quantum 
computation" O 0] employing non-Abelian statistics of some of the wave functions proposed 
for a partially filled LLi. The main idea is that the quantum information may be encoded in 
(topologically) different quantum states corresponding to the same spatial configuration of the 
non-Abelian quasiparticles (QPs) of a given underlying IQL. Transition between such different 
states would only be possible by a global transformation of the QP braiding, thereby making 
quantum information inherently protected from decoherence caused by any local processes (e.g., 
coupling to phonons or atomic spins). The best known candidate for a non-Abelian wave function 
is the "pfaffian" state proposed by Moore and Read E] and believed to describe the FQH 
state in a half-filled LLi. Other wave functions with different complexities of braiding statistics 
have also been proposed El [TU] . The convincing demonstration of non-Abelian statistics in 
a real physical system has clearly risen to the challenge of greatest importance. 

The reason to search for the signatures of non-Abelian statistics in LLi is that, on the one 
hand, the partially filled LLo are successfully described by the composite fermion (CF) theory 
predicting fractional but Abelian QPs and, on the other hand, that higher LLs favor ordered 



electron phases. Crucial recent experiments in LLi include confirmation [12] of the anticipated 
QP charge of e/4 for the half- filled state at v e = 5/2 and careful measurements |13|. [T3] of the 
minute excitation gaps. In theory, some of the most recent advances are related to the role of 
layer width [15] and LL mixing [161 117 j in real systems, and the nature of the QPs |18t I19j. 

It is quite remarkable that, despite intensive studies, connection of the FQH states in LLi 
[y e = 5/2, 7/3, 12/5, or 11/5) to the few proposed wave functions remains yet to be conclusively 
established. In fact, in some cases it seems only tentatively assumed for the lack of other 
candidates. This is an urgent problem, as the anticipation of non-Abelian statistics in nature 
is largely driven by the connection of some of these wave functions to the particular conformal 
field theories. A wealth of IQLs found in various systems (electrons or CFs at different fillings of 
different LLs, in layers of varied width w) also invites a more general question of possible IQLs 
with arbitrary LL filling v and interaction V. 

This paper is an extension of our brief communication [20] . We report here on the use of large 
scale exact diagonalization in the search of IQLs with non-Abelian statistics. We demonstrate 
that non-Abelian IQLs in LLi can only emerge in the narrow, particle-hole symmetric range of 
filling factors 7/3 < u e < 8/3. In this range, the known non-Abelian candidates at v e = 5/2 and 
12/5 are closely examined (including such previously neglected realistic effects as finite layer 
width and LL mixing) and found to have favorable correlation energies. Outside of this range, 
the LLi hosts the family of Abelian ground states of noninteracting CFs each carrying four 
magnetic flux quanta, repeating the states known from the lowest LL. The borderline v = 1/3 
ground state in LLi is adiabatically connected to the Laughlin state of LLo, but it has a smaller 
gap and distinct short-range correlations. 

2. Model 

We consider the systems of N spin-polarized fermions (electrons) on a Haldane sphere with unit 
radius and the magnetic monopole of strength 2(^(/ic/c) inside [21j . In this geometry, LL n is 
a shell of angular momentum I = Q + n, and different iV-body wave functions at the same v 
are (unlike on a torus) conveniently distinguished by a 'shift' 7 between the LL degeneracy and 
v v N (i.e., 21 = v~ 1 N — j). In contrast to the previous exact diagonalization studies we begin by 
searching the universality classes (u, 7) of the gapped ground states with arbitrary interactions 
V rather than confining ourselves to the particular physical systems defined by n, uu, etc. 

It is a trivial fact that the many-body dynamics in a degenerate LL is completely determined 
by an interaction pseudopotential, defined [21] as the dependence of the pair energy V on the 
relative angular momentum m = 1, 3, .... Less obviously, V m induces particular correlations 
only through its deviation from a reference "harmonic pseudopotential" given by a straight line 
over the consecutive m's [22 j . Hence, the low-energy spectra of V m can be accurately reproduced 
by a suitable effective pseudopotential U m with only a few appropriate coefficients. 

In our calculations we have used U = [U\, U3, U5]. Higher-order terms have been ignored, as 
they are essentially irrelevant for the dynamics in a liquid phase with short-range correlations. 
On the other hand, going beyond the (earlier used) U\ and U3 was needed for an improved 
description of the two lowest (electron or CF) LLs known to host IQL states. At the same time, 
it still allowed for useful graphical representation of the ground state properties in an effectively 
two-dimensional space of the (normalized, ^ m U m = 1) parameters U m . 

3. Maps of the gap for arbitrary interaction 

Searching for the series of gapped ground states with particular filling factor v and shift 7, we 
looked at various finite systems (iV, 2£). A few maps of the 'neutral' energy gap A for v = 1/2 
and 7 = 3 (i.e., 21 = 2N — 3, as appropriate for the Moore-Read pfaffian state) are shown in 
Fig. [TJ The gap A is defined as the energy difference from the ground state to the first excited 
state in the same spectrum as long as the ground state happens to be nondegenerate (i.e., has 




Figure 1. Ternary contour plots of the neutral excitation gap A for N = 10, 12, 14, 16 fermions 
in a half-filled Landau level (y = 1/2) with shell angular momentum £ corresponding to the shift 
7 = 3 of the Moore-Read "pfaffian" wave function, interacting by model pseudopotential U m . In 
each plot, three corners of the big triangles correspond to U m = 5 m ^ with = 1, 3, 5 marked in 
each corner. Points relevant for the actual interactions in different electron or CF Landau levels 
are marked on a triangle in frame (e). Different candidate incompressible states are indicated. 



zero total angular momentum, L = 0); otherwise A is set to zero. The IQL candidates are 
the islands of positive A repeating regularly in the same area of the map for different values of 
N. Their location on the map must be compared with the actual pseudopotentials in different 
LLs, as indicated in frame (e). In particular, (A) and (B) mark the appropriate positions of 
the electron pseudopotentials in LLo and LLi (the former dominated by V\\ the latter roughly 
linear between m = 1 and 5), and (X), (Y), and (Z) mark the same for the CFs in their effective 
shells LLq, LLJ, and LL^. 

The appearance of a significant gap around the point (B) in all maps in Fig. [T] (i.e., for 
each N) confirms quite definitively the earlier expectation that the Moore-Read ground state 
forms for a class of pseudopotentials close to that of LLi. It also demonstrates that its accuracy 
depends sensitively on the fine-tuning of the leading V^'s, achieved (for example) by adjusting 
the layer width w [15] . Remarkably, the maps in Fig. Q] also preclude the Moore-Read state at 
the half-filling of other LLs (e.g., of the second CF LL, called LL*, characterized by a dominant 
repulsion at m = 3 [23]). Clearly, the microscopic origin of the FQH state observed at v e = 3/8 
[24] (in the CF picture, corresponding to v = 1/2 in LLJ) must be different. 

Analogous maps of A for the Laughlin filling fractions v = 1/3 and 1/5 are shown in Fig. [2j 
In (a), the universality class of the Laughlin wave function {y = 1/3 and 7 = 3) is shown to 
cover a large part of the map, including (A) and (X), and possibly also (B) and (Z). Point 
(B), the most interesting for the present analysis, falls just inside the island of positive gap A, 
suggesting connection of the v e = 7/3 FQH state to the Laughlin v = 1/3 liquid in LLi. Point 
(X) is relevant for the Laughlin v = 1/3 state of the CF vacancies in LLq, i.e., to the robust 



Figure 2. Similar to Fig. [TJ but for Laughlin 
filling fractions v = 1/3 (two candidate states 
with 7 = 3 and 7) and v = 1/5. 



Figure 3. Similar to Fig. [TJ but for the v - 

2/5 states with 7 = 4 (Jain), 2 (Bonderson 
Slingerland) , and —2 (Read-Rezayi) . 



Jain state at v e = 2/7. Remarkably, the Laughlin ground state does not occur around point 
(Y) corresponding to the v = 1/3 filling of LL*. The nature of the rather fragile FQH state 
observed at the corresponding fraction v e = 4/11 |24j must therefore be different. In (b), the 
v = 1/3 paired (non-Laughlin) state with 7 = 7, proposed earlier for both LLi [22J and LL^ [25], 
is tested. The gap around (Y) is quite suggestive that it may indeed describe the FQH state at 
v e = 4/11 [23]. On the other hand, its relevance for LLi seems doubtful. In (c), any positive 
U = [U\, U3, 0] yields an exact Laughlin state at v = 1/5. Our map of A confirms that it is true 
description of the FQH states in both LLo and LLi [26]. On the other hand, its relevance to the 
FQH effect observed in LLJ (i.e., at v e = 6/17) [23] is rather unlikely. 

Finally, the maps of A for the v = 2/5 filling have been shown in Fig. [3j In (a), the Jain 
series of noninteracting CF states with 7 = 4 correctly represents the ground state in LLo (and 
LLj$). Point (B) corresponding to LLi lies just outside the borders of the island of A > 0, 
leaving the question of relevance of the non-interacting CF picture at v e = 12/5 open. Two 
other candidate ground states at v = 2/5 are the parafermion state with 7 = — 2 [7] and a 
more recently proposed state with 7 = 2 [9j. Especially for the latter state, frame (b) appears 
suggestive of a gap emerging around (B). Clearly, the competition between these three candidate 
ground states in LLi is not convincingly resolved based on the maps of A alone. On the other 
hand, the identification of the true ground state is crucial, because the candidates with 7 = ±2 
are both non-Abelian, in contrast to the 7 = 4 Jain state. More careful analysis will follow in 
subsequent sections. 

4. Maps of the amplitudes for arbitrary interaction 

Maps similar to those in Figs. [TH21can also be used to show dependence of other parameters of the 
spectrum on the form of interaction. For example, in Fig. [3] we plot maps of the leading pair and 
triplet Haldane amplitudes for the ground states of [y, 7) = (1/2, 3) and (1/3, 3). The amplitudes 
are defined [2Z] as the fractions of the pairs or triplets with a given relative angular momentum m 
(for the pairs, m = 1, 3, 5, . . . is a measure or an average square distance; for the triplets, m = 3, 

5, 6, . . . measures an average area). These (discrete) correlation functions are particularly useful 
in identifying Laughlin and Moore-Read states, as they both are unique zero-energy states of 
simple repulsions: the former at the minimum pair angular momentum m = 1, the latter at 
the minimum triplet angular momentum m = 3 (hence, their corresponding amplitudes vanish 
exactly). The emergence and location of the islands of essentially zero amplitude in Fig. |3]^c) and 
(d) provides additional and quite decisive support for the Laughlin v = 1/3 and Moore-Read 
v = 1/2 ground states (universality classes) in LLo and LLi, respetively. On the other hand, 
other questions, such as of a Laughlin state in LLi, remain open until the competition with 
other possible states (with the same v but different 7) can be resolved. 
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Figure 4. Ternary contour plots (similar to Fig. [T]) of the pair and triplet amplitudes (labeled 
by the relative angular momentum m) for the universality classes of Laughlin and Moore-Read 
wave functions (filling factors v = 1/3 and 1/2, respectively; shift 7 = 3 in both cases). 




Figure 5. Symmetric second-order differences 5 2 E of the ground-state energy per particle E for 
N = 10 and 12 electrons in the lowest and second Landau level (LLo and LLi), as a function of 
shell angular momentum £, for layer widths w = and 3A. Candidate incompressible states are 
marked as where $ = P, P, L, J, B denote Moore-Read pfaffian, anti-pfaffian, Laughlin, 
Jain, and Bonderson-Slingerland states. Squared overlaps \\\ 2 between the states repeating in 
both LLs at v < 1/3 are indicated. A is the magnetic length. 



5. Analysis of correlation energies in LLo an d LLi 

Guided by the maps of gaps and amplitudes we now move our focus to the FQH states in a 
partially filled LLi. In Fig. [5] we seek confirmation of the IQL candidates in the downward cusps 
of the dependence of the ground-state correlation energy per particle E on the LL degeneracy 
(for a fixed number of electrons N). The correlation energy E is calculated from the total 
Coulomb energy £ of N electrons in a LL shell with a given £ by adding the energy of attraction 
to the uniform charge-compensating background and dividing by N. In the calculation for finite 
width w of a quasi-2D electron layer, the Coulomb matrix elements were computed assuming 
infinite- well confinement in the perpendicular direction, i.e., for the charge-density profile of 
g(z) = (2/w)cos 2 (ttz/w). The cusps in E are most pronounced in the plots of a symmetric 




Figure 6. Size extrapolation of the ground-state correlation energies per particle E calculated 
for N electrons in the second Landau level (LLi), for the layer widths w = (a) and 3A (b). 
Competing series of candidate incompressible ground states with filling factors v = 1/3 and 2/5 
are distinguished by their shifts 7. 



difference 5 2 E2e = E^t-x + -EW+i — 2E2e- For an IQL, its (positive) value gives the QP gap A 
(energy needed to create a pair of noninteracting QPs of total charge zero) times the number of 
QPs created per flux quantum. 

Evidently, Fig. [5] complements the maps of Fig. [T0 in the identification of the universality 
classes of particular IQLs. For example, it shows peaks in 5 2 E2e which signal the Laughlin 
v = 1/3 and Moore-Read v = 1/2 IQLs in LLi. But Fig. [5] also does more, by revealing the 
following connection between quantum statistics and the filling factor in a partially filled LLi. 
At v < 1/3 the same Laughlin and Jain IQLs of filled shells of noninteracting CFs with four 
flux quanta occur in LLi and LLo- The corresponding states in LLi and LLo have high overlaps 
(calculated by replacing the electron positions by the guiding centers) and similar QP gaps A. 
This similarity, earlier pointed out in Ref. [20], is caused by a sufficiently high pseudopotential 
coefficient V\ (forcing the avoidance of the m = 1 pair state) and a similar behavior of the 
pseudopotential al long range, V m >s, in the two lowest LLs. Importantly, this similarity validates 
the noninteracting CF model [TT] with four flux quanta attached to each electron in LLi (in 
addition to LLo where its accuracy is well-known). 

In contrast, at v > 1/3 the Jain sequence of noninteracting CF states in LLo is replaced in 
LLi by a different set of IQLs, including the Moore-Read pfaffian at v = 1/2, the anti-pfaffian 
(pfaffian's particle-hole conjugate) at the same u = 1/2 but with a different shift 7 = — 1, and a 
Bonderson-Slingerland v = 2/5 state with 7 = 2. The breakdown of the noninteracting two-flux 
CF model in LLi opens, exclusively at 1/3 < v < 2/3, possibility for other IQLs, including 
several suggested more exotic states with various non-Abelian QP statistics. The borderline 
v = 1/3 state, separating the Abelian from (possibly) non-Abelian states, has a moderate 
overlap with the Laughlin state of LLo, despite falling into the same class of 7 = 3. 

Having narrowed the search for non-Abelian IQLs in LLi to the range of 1/3 < v < 2/3, let 
us now try to resolve more decisively the competition between different candidate wave functions 
at filling factors v = 1/3 and v = 2/5 (at v = 1/2, the Moore-Read pfaffian appears to have 
no competition). In Fig. [6] we attempt to extrapolate to an infinite system size (to iV^ 1 — » 0, 
i.e., to the planar geometry) the ground-state correlation energies per particle E. To improve 
convergence, the plotted energies E have been rescaled by \/2Qv/N so as to ensure equal units 
e 2 /A for each 7; here A = y/hc/eB is the magnetic length). In the extrapolation for u = 2/5 we 
only used the open squares, discarding one state aliased with the anti-pfaffian (a finite system 
defined by N and 21 may sometimes represent states with different combinations of v and 7) 



and one apparently suffering from the small size. When choosing the IQL ground state we also 
checked if a given candidate series (i/, 7) consistently has a nondegenerate (i.e., one with L = 0) 
ground state. In those cases when the particular (N, 21) ground state had L / 0, we crossed 
out the corresponding data point in Fig. [6] to mark that the series containing it is unlikely to 
describe an IQL in the thermodynamic limit. 

For v = 1/3 Fig. [6] eliminates the (paired) 7 = 7 series as a viable candidate in LLi (it remains 
plausible in LL| which, however, will not be further discussed here), leaving the (universality 
class of the) Laughlin state as the only acceptable option. For v = 2/5, the Jain series appears 
to extrapolate to a competitive energy, but it no longer has L = for the larger N, and hence 
should probably be discarded. The remaining two candidates both consistently have L = and 
both extrapolate to nearly the same E for large iV (the difference in favor of the Read-Rezayi 
state with 7 = —2 is comparable to the error of extrapolation). In realistic conditions this near 
degeneracy is likely to be removed by their different susceptibilities to the LL mixing. To explore 
this idea, we estimated the appropriate energy correction dE by including in the diagonalization 
additional states involving a single cyclotron excitation [16J. Specifically, in addition to the 
configurations with completely filled LLo (with both spin-f and spin-j), ./V polarized electrons 
in LLi, and empty higher LLs, we have also included configurations with one electron promoted 
from either LLo (with either spin) or from LLi to the next higher LL. It should be kept in mind 
that inclusion of only a single cyclotron excitation is not a rigorous treatment of the LL mixing. 
However, it is expected to reveal a possible difference in the susceptibility of the competing states 
to this process. In the calculation, we have assumed the Coulomb-to-cyclotron energy ratio of 
(3 = (e 2 / \)/(hu) c ) = 1.56, corresponding to B = 2.6 T. Due to a large size of the Hilbert space, 
the values of dE could only be calculated for N < 10. However, we found that the corrections 
dE are far less size-dependent than the base energies E. It is therefore justified to apply the 
N = 10 estimates of dE to the values of E extrapolated from N < 16. 

The results are presented in Tab. [TJ Clearly, in addition to having the lowest E, the 7 = —2 
parafermion state also has the largest \dE\, and therefore it is quite convincingly predicted 
to define the universality class of the v = 2/5 ground state in LLi. This is of considerable 
importance because this state is the only known candidate IQL [4] whose braiding rules are 
sufficiently complex to allow quantum computation (unlike, e.g., the Moore-Read state). 

Let us now look at the pair correlation functions of the established IQL ground states at 
v = 1/3 and 2/5 in LLi. The results calculated for the largest available systems (N = 14 or 16, 



Table 1. Extrapolated correlation energies per particle E (in the units of e 2 /A) in several series 
of candidate ground states (distinguished by shifts 7) at filling factors v = 1/3 and 2/5 in LLi, 
for layer widths w = and 3A, calculated without LL mixing. Also, corrections dE due to LL 
mixing (at B = 2.6 T), estimated for N = 10 (except for N = 8 for v = 2/5 and 7 = 4). 



0, 7) 


(1/3,3) 


(1/3,7) 


(2/5,-2) 


(2/5,2) 


(2/5,4) 


w = 0: 
E 
dE 
E + dE 


-0.3285 
-0.0268 
-0.3553 


-0.3249 
-0.0218 
-0.3467 


-0.3439 
-0.0294 
-0.3733 


-0.3406 
-0.0239 
-0.3645 


-0.3427 
-0.0227 
-0.3654 


w = 3A: 
E 
dE 
E + dE 


-0.2611 
-0.0087 
-0.2698 


-0.2589 
-0.0077 
-0.2666 


-0.2657 
-0.0095 
-0.2752 


-0.2642 
-0.0085 
-0.2727 


-0.2659 
-0.0071 
-0.2730 




Figure 7. Pair-correlation functions of the candidate incompressible ground states at v = 1/3 
(a) and 2/5 (b) in the second Landau level (LLi). Coulomb ground states are compared with 
the exact Laughlin and approximate Jain states in LLo and LLi (constructed as the ground 
states of a model interaction pseudopotential with pair repulsion only at m = 1). 




Figure 8. Renormalized pair amplitudes T of the same candidate states in LLi as in Fig. [JJ 



depending on a given state, as is clear from Fig. [6]) are shown in Fig. [JJ At v = 1/3, despite 
being adiabalitally connected to the Laughlin state (i.e., having the same 7 = 3), the Coulomb 
ground state has significantly different short-range correlations. Specifically, the zigzag in g{r) 
atl<r/A<2 present in the exact Coulomb state in LLi, caused by the particular form of the 
single particle wave functions and showing also in g(r) for the full LLi, gives way to a plateau 
in the case of the actual Coulomb ground state. (We should also clarify here that the Laughlin 
state in any LL is defined in a standard way as the zero energy state of the m = 1 pair repulsion, 
and its pair correlation function is drawn using the appropriate single-particle wave functions, 
e.g., of LLo ° r LLi.) This agrees with only moderate squared overlaps of the v = 1/3 ground 
state in LLi with the Laughlin wave function, pointed out ourlier: 0.292 (0.501), 0.253 (0.510), 
0.333 (0.549) for N = 10, 12, 14, for w = (3A). The magneto-roton band is also absent in 
the spectra of LLi, and the low-energy states resembling Laughlin QEs and QHs are found at 
21 = 3N — 3^1, but they are not generally the lowest states in their spectra. A similar difference 
(removal of the characteristic zigzag) between the Jain state (approximated as the ground state 
of the m = 1 pair repulsion for 7 = 4) and the competing Coulomb ground states with 7 = —2 
and 2 is also found at v = 2/5 in LLi. At both filling factors dependence on w is insignificant. 

Difference in short-range correlations between Laughlin and Jain states of LLi on one hand, 
and the v = 1/3 and 2/5 Coulomb ground states of LLi on the other, was also earlier apparent 



in the leading amplitudes in Fig. Efd,e). It is evident in the comparison of renormalized 
pair amplitudes T m calculated for slightly larger systems (same as in Fig. UJ), which have 
been shown in Fig. [8j Here, the coefficients T m have been converted from the usual Haldane 
amplitudes on a sphere G m to their planar counterparts, as described in Ref. [27]. Specifically, 
r m = 1 + (G m - G^ ll )/G[ ul \ where G^ 11 = (4£ + l- 2m) /I/ (21 + 1) describes a full LL. 



6. Conclusion 

We have carried out extensive exact diagonalization studies (including finite layer width and LL 
mixing) of the fractional quantum Hall states in a partially filled second Landau level (LLi), 
searching for non-Abelian incompressible quantum liquids (IQLs). We have found the range of 
filling factors 1/3 < v < 2/3 in LLi in which the emergence of non-Abelian statistics is possible. 
Inside this range, we have demonstrated that the spin-polarized ground states at v = 1/2 and 
2/5 are described by the non-Abelian "pfaffian" and "parafermion" wave functions. Outside of 
it, the Jain states of noninteracting CFs with four attached magnetic flux quanta repeat in both 
lowest LLs, thus precluding more exotic phases. The borderline u = 1/3 state is adiabatically 
connected to the Laughlin liquid but has a smaller gap and distinct short-range correlations. 
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